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Abstract 

Porcine sapelovirus (PSV), a species of the genus Sapelovirus within the family Picomaviridae, is associated with diarrhea, 
pneumonia, severe neurological disorders, and reproductive failure in pigs. However, the structural features of the complete 
PSV genome remain largely unknown. To analyze the structural features of PSV genomes, the full-length nucleotide 
sequences of three Korean PSV strains were determined and analyzed using bioinformatic techniques in comparison with 
other known PSV strains. The Korean PSV genomes ranged from 7,542 to 7,566 nucleotides excluding the 3' poly(A) tail, and 
showed the typical picornavirus genome organization; 5'untranslated region (UTR)-L-VP4-VP2-VP3-VP1-2A-2B-2C-3A-3B-3C- 
3D-3'UTR. Three distinct c/s-active RNA elements, the internal ribosome entry site (IRES) in the 5'UTR, a c/s-replication 
element {CRE) in the 2C coding region and 3'UTR were identified and their structures were predicted. Interestingly, the 
structural features of the CRE and 3'UTR were different between PSV strains. The availability of these first complete genome 
sequences for PSV strains will facilitate future investigations of the molecular pathogenesis and evolutionary characteristics 
of PSV. 
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Introduction 

Picornaviruses are a family of positive-sense single stranded 
RNA viruses within the order Picornavirales [1]. They can cause 
intestinal, respiratory, neurological, cardiac, hepatic, mucocuta- 
neous, and systemic diseases of varying severity in humans and 
animals [2]. Although different picornaviruses show various 
degrees of relatedness, all picornaviruses share a similar genomic 
organization, which consists of a covalently linked 5' terminal 
protein called VPg (Viral Protein genome-linked), a 5' untrans- 
lated region (UTR), a large open reading frame (ORF), a 3' UTR 
and a poly(A) tail of variable length [2,3]. The genomic RNA of 
picornaviruses harbor several distinct CM-active RNA elements 
which are required for viral RNA replication; an internal ribosome 
entry site (IRES) in the 5'UTR, a cm -replication element (CRE) 
within the ORF [3-6] or the 5' UTR [7], the 3' UTR, and the 3' 
poly(A) tail [3]. Currently, five different types of IRES element [8] 
that direct cap-independent translation initiation on the viral RNA 
to produce the polyprotein have been identified from the primary 
sequence, secondary structure, location of the initiation codon and 
activity in different cell types [9,10]. In most picornaviruses, the 
polyprotein encoded by the ORF is cleaved into four structural 
viral particle proteins (VP4, VP2, VP3 and VP1) and seven non- 



structural proteins (2A, 2B, 2C, 3A, 3B, 3C and 3D). In addition, 
the members of the genera Cardiovirus, Aphthovirus, Erbovirus, 
Kobuvirus, Teschovirus, Senecavirus and Sapelovirus possess a 
leader (L) protein at the N-terminus of polyprotein [11]. 

Although simian type 2 picornaviruses (SV-2-like viruses) and 
porcine enterovirus 8 (PEV-8) were once classified in the genus 
Enterovirus, SV-2-like viruses and PEV-8 along with duck 
picornavirus TW90A have an L protein at the N-terminus of 
the polyprotein that is lacking in the enteroviruses [12—14]. 
Moreover, those viruses contain distinct 2A proteins from those of 
the Enterovirus genus, and a highly divergent 5'UTR with a type 
IV IRES [12,14-16]. Due to these particular genetic properties, 
these simian, avian, and porcine picornaviruses are now assigned 
as members of a new picornavirus genus, Sapelovirus [12— 
15,17,18]. 

PSV infections have been associated with a wide spectrum of 
symptoms ranging from asymptomatic infection to clinical signs 
including diarrhea, pneumonia, polioencephalomyelitis, and 
reproductive disorders [19-20]. Although PSVs can be important 
pathogens because of their wide distribution and high prevalence 
[21-25], the near-complete genomic sequences of only three PSV 
strains have been reported previously; one from the U.K. (VI 3 
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strain) and two from China (csh and YC2011 strains). This 
prompted us to characterize the full-length genetic properties of 
Korean PSV strains in comparison with those of other known PSV 
strains. 

During an epidemiological study on PSV infections in the fecal 
samples of piglets with diarrhea in South Korea, three PSV strains 
were isolated. These Korean PSV strains were characterized using 
an immunofluorescence assay (IFA) with a monoclonal antibody 
specific for a PSV capsid protein, RT-PCR assay with primer pair 
specific for the PSV VP 1 coding region and transmission electron 
microscopy. Furthermore, bioinformatic techniques were em- 
ployed to analyze the complete viral genomes of the three newly 
isolated strains in comparison with the other known PSV strains. 

Materials and Methods 

Origin of virus strains 

The diarrheic fecal samples, collected from piglets in 45 
different herds in South Korea during 2004—2007, were screened 
for PSV infections using RT-PCR and nested-PCR assays with 
primer pairs specific for the PSV VP1 coding region [26]. Among 
the PSV-positive fecal samples, three fecal samples which were 
strongly positive for PSV by RT-PCR were used to isolate PSVs 
using a porcine kidney cell line [22], LLC-PK. In brief, the RT- 
PCR positive fecal samples were diluted 10-times with 0.01 M 
phosphate-buffered saline (PBS, pH 7), vortexed for 30 s and 
centrifuged at 1200xg for 20 min. The supernatants were filtered 
through 0.2-um syringe filters. Filtered supernatants were serially 
diluted 10-times with Eagle's minimal essential medium (EMEM) 
containing 1% antibiotics (Penicillin, Streptomycin, and Ampho- 
tericin B) and 1% NaHC0 3 , and used to infect cells in 6-well 
plates. The suspensions were absorbed for 1 h with occasional 
rocking, and EMEM containing 1 % antibiotics and 1 % NaHCO :i 
was added. The cultures were incubated for 3 to 4 days at 37°C in 
a 5% C0 2 atmosphere and examined daily for cytopathic effects 
(CPE). Isolated PSVs were cloned by triple plaque purification. 
The PSV strains (KS04105, KS05151 and KS055217) were 
passaged eight times in LLC-PK cells, including isolation, 
adaptation, and triple plaque purification. The isolated viruses 
were confirmed as PSVs by an IFA, RT-PCR and transmission 
electron microscopy (TEM) assays, as described below. 

Transmission electron microscopy (TEM) 

LLC-PK cells infected with each of the above strains and 
showing over 70% CPE were frozen and thawed thrice, and 
centrifuged at 2,000 xg for 30 min. To obtain purified virus, each 
supernatant was ultra-centrifuged at 200,000 xg for 5 h at 4°C in a 
S58A-0015 rotor (Hitachi, Tokyo, Japan). The resulting pellets 
were resuspended in 40 ul of water and mixed with an equal 
volume of 2% (w/v) sodium phosphotungstic acid at pH 7.0. The 
samples were placed onto a formvar grid (Electron Microscopy 
Sciences, Hatfield, USA) for 5 min, and then excess liquid was 
removed by filter paper. The samples were examined using a High 
Resolution Transmission Electron Microscope (Hitachi) for the 
determination of purity of virus stock at Gwangju Center of Korea 
Basic Science Institute. 

Immunofluorescence assay (IFA) 

To characterize the PSV strains, the IFA was performed with a 
153/5B5 (IgG2a) monoclonal antibody specific for the PSV capsid 
protein (kindly provided by Dr. M Dauber, Friedrich-Loeffler- 
Institut, Greifswald-Insel Riems, Germany). Briefly, LLC-PK cells 
were infected with each strain, incubated for 18 h as above, fixed 
in 80% acetone for 5 min at 4°C, and allowed to air dry. Slides 



were washed thrice with PBS (pH 7.2), and incubated overnight at 
4°C using a 1:40 dilution of monoclonal antibody specific for PSV 
capsid protein diluted in PBS (pH 7.2). Slides were washed thrice 
with PBS (pH 7.2), and incubated with FITC-conjugated goat 
anti-mouse IgG antibody (Santa Cruz biotechnology, Santa Cruz, 
USA) diluted 1:100 in PBS (pH 7.2) for 1 h at room temperature. 
After washing twice with PBS (pH 7.2), slides were stained with 
4',6-diamidino-2-phenylindole (DAPI) (Invitrogen, Lohne, Ger- 
many), and examined using a LSM confocal scanning microscope 
(Carl Zeiss, Jena, Germany). 

RNA extraction and RT-PCR 

Total RNA was extracted from the lysates of LLC-PK cells 
infected with each strain using the AccuPrep Viral RNA 
extraction kit (Bioneer, Daejeon, Korea) according to the 
manufacturer's instructions. To detect and amplify PSV RNA, 
RT-PCR with a primer pair specific for the PSV VP1 coding 
sequence (Table SI in File SI) was performed. To characterize the 
complete full-length genome sequences of each strain, ten primer 
sets (Table SI in File SI) were designed to amplify the complete 
ORF sequences of each PSV strain based on the published 
genomic sequences of the PSV-V13 (GenBank ID: NC_003987), 
csh (GenBank ID: HQ875059) and YC201 1 strains (GenBank ID: 
JX286666). Standard one-step RT-PCR assays were performed as 
previously described [25]. 

5' and 3' cDNA syntheses 

cDNA of each strain was synthesized by the SMARTer Rapid 
Amplification of cDNA Ends (RACE) cDNA amplification kit 
(Clontech, Mountain View, USA) according to the manufacturer's 
instructions. For generating 3' RACE-ready cDNA, 3.75 ul of the 
poly(A) tailed RNA and 1 ul of 3'-cDNA Synthesis (CDS) Primer 
A were mixed and heated to 72°C for 3 min, followed by cooling 
to 42°C for 2 min using a thermo cycler. For generating 5' 
RACE-ready cDNA, 3.75 ul of total RNA was mixed with 1 ul of 
5'-CDS primer A, incubated at 72°C for 3 min, and cooled at 
22°C for 2 min. The denatured RNA for each 3' and 5' cDNA 
generation was mixed with a reaction mixture composed of 2 u.1 
5 x First-Strand Buffer, 1 ul dithiothreitol (DTT) (20 mM), 1 ul 
dNTP mix (10 mM), 0.25 ul RNase inhibitor (40 U/ul), and 1 ul 
SMARTScriber Reverse Transcriptase (100 U). Samples were 
incubated at 42°C for 90 min and heated at 70°C for 10 min. The 
synthesized cDNAs were diluted with 7 ul of Tricine-EDTA buffer 
and used for RACE PCR. 

RACE PCR, cloning and sequencing 

For the generation of 3' and 5' RACE PCR reactions, 
Advantage 2 Polymerase Mix (Clontech) was used; 5 ul of 3'/5'- 
RACE-Ready cDNA, 32 ul of PCR-grade water, 5 ul of 
lOxAdvantage 2 PCR Buffer, 1 ul dNTP Mix (10 mM), 5 ul 
Universal Primer Mix (10 X), 1 Jul of 50 pmol/ml gene-specific 
primer (GSP) for 3' and 5' RACE (Table SI in File SI), and 1 ul 
of 50 xAdvantage 2 Polymerase Mix. The reaction was performed 
with the following thermal cycling program: 5 cycles of 94°C or 
94.5°C for 30 sec and 72°C for 2 min or 3 min; 5 three-step cycles 
of 94°C or 94.5°C for 30 sec, 65°C (applied melting temperature 
of GSP) for 30 sec, and 72°C for 2 min or 3 min; 25 three-step 
cycles of 94°C or 94.5°C for 30 sec, 60°C (applied lowered Tm 
values of GSP by 3°C to 5°C) for 30 sec, and 72°C for 2 min or 
3 min. 

The RACE PCR products were separated on a 1.5% (m/v) 
agarose gel; the bands were excised and purified using a Purelink 
Quick Gel Extraction kit (Invitrogen). The products were ligated 
into TA Vector Systems (Enzynomics, Daejeon, Korea) and 
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introduced into DH5a competent cells. Individual colonies were 
grown and plasmid was purified using Hybrid-Q_ Plasmid 
(GeneAll, Seoul, Korea). Sequencing was performed using an 
ABI System 3700 automated DNA sequencer (Applied Biosystems, 
Foster, USA). 

Full-length genomic characterization and secondary 
structure prediction 

To generate complete nucleotide sequences for each PSV strain, 
both 5 ' and 3 ' end sequences of each strain were assembled with 
the nucleotide sequences of the internal ORF sequences. The 
complete full-length genomic and individual protein coding 
sequences of three PSV strains were compared with those of the 
other known PSV strains (Table S2 in File SI) using the DNA 
Basic module (DNAsis MAX, Alameda, USA). Phylogenetic 
analyses based on nucleotide and amino acid alignments were 
performed using the neighbor-joining method with 1000 bootstrap 
replicates and UPGMA Molecular Evolutionary Genetics Analysis 
(MEGA version 5.2) employing pair-wise distance comparisons 
[27]. Sequence identity calculations for the three PSV strains with 
those of the other known PSV strains were performed using the 
homology and distance matrices method of DNAMAN version 6.0 
program (Lynnon, Vaudreuil, Canada). Secondary structure 
elements in the PSV genomes were modeled using GLC Main 
Workbench version 6.8.2 program (CLC bio, Katrinebjerg, 
Denmark). 

Ethics statement 

No specific approval was needed since the fecal samples were 
voluntarily submitted by the farms for pathogen screening in our 
laboratory. No other specific permits were required for the 
described field studies. The locations where we sampled are not 
protected in any way. The field studies did not involve endangered 
or protected animal species. Before beginning work on the study, 
we contacted the farm owners and obtained their permission. 

Results 

Virus isolation and identification 

Three Korean PSV strains were isolated from separate diarrhea 
fecal samples originating from three different farms and plaque 
purified. LLC-PK cells infected with each strain after 8 passages in 
LLC-PK cells showed CPE at day 1 post-inoculation characterized 
by shrinkage, rounding and detachment of cells (Figure 1A, IB), 
and displayed PSV-specific cytoplasmic fluorescence in the 
indirect IFA using a monoclonal antibody against PSV capsid 
protein (Figure 1 A). RT-PCR assays with a primer pair specific for 
the partial PSV VP1 coding region amplified a 636 bp amplicon 
from LLC-PC cells infected with each strain (Figure SIC). By 
transmission electron microscopy, negatively stained purified virus 
particles of each strain appeared spherical with a diameter of 
approximately 30 nm (Figure IB). No other virus-like particles 
were observed. These results identified the isolated viruses as 
PSVs. 

Genome organization 

The complete nucleotide sequences of the whole genome of the 
three Korean PSV strains were obtained and compared to the 
previously determined PSV sequences. In picornavirus RNAs, the 
two 5' terminal UU residues are derived from the uridylylation of 
VPg to make VPg-pU-pU [28]. In the previously described PSV 
sequences, however, the two 5' terminal residues were AC for the 
Chinese csh and English VI 3 strains, and UA for the Chinese 
YC201 1 strain, suggestive of incomplete sequences. To obtain the 



correct 5' terminal start residues, cDNA synthesis and then 5' 
RACE PCR with 5' RACE primer (Table SI in File SI) were 
performed. Using this approach, the 5' terminal residues were UU 
and the 5' UTR length of three Korean strains was 25 nucleotides 
longer than that of Chinese YC2011 strain (Table S3 in File SI) 
[29]. In order to confirm this result, 5 'RACE PCR was performed 
with another 5'RACE PCR primer (Table SI in File SI) and the 
same 5' terminal nucleotide residues, UU, were also observed. 

The length of the complete genomes of Korean PSV strains, 
excluding the poly(A) tail, was from 7,542 (KS04105 and 
KS055217) to 7,566 nucleotides (KS05151) (Table S3 in File 
SI). These sequences contained a single large ORF whose lengths 
were 6,966 nucleotides (strain V13; 2,322 amino acid polyprotein 
precursor) to 6,993 nucleotides (strains csh, KS05151, YC2011; 
2,331 amino acid polyprotein precursor) (Table S3 in File SI). The 
predicted protease cleavage sites of the polyproteins, as determined 
from alignments with other picornaviruses are shown in Table 1. 
The polyprotein was predicted to be cleaved and processed into 
twelve mature peptides: L-VP4-VP2-VP3-VP1-2A-2B-2C-3A-3B- 
3C Pr "-3D Po1 (Figure 2). The ORF sequence in the three Korean 
PSV strains was flanked by 5' UTR which was 491 nucleotides 
long and by a 3' UTR which was 82 nucleotides long (Table S3 in 
File SI). 

Molecular and phylogenetic analyses 

The complete genome sequence, excluding the poly(A) tail, and 
the polyprotein sequences of three Korean strains were compared 
with those of other known PSVs and representative picornavirus 
strains available in the GenBank database. The Korean PSV 
strains showed high nucleotide (84.7%-94.0%) and deduced 
amino acid (92.9%-98.3%) identities with the other PSV strains 
(Table 2), but showed relatively low nucleotide and polyprotein 
sequence identities with the avian and simian sapelovirus strains 
(Table 2). 

Each of the major functional units in the genome, including the 
5' and ?>' UTRs, the capsid coding region (PI) and the regions 
encoding the non-structural proteins (P2 and P3) of PSVs, were 
phylogenetically analyzed (Figure 3). Representative simian and 
avian sapelovirus reference strains were included in each of the 
trees. The 5' UTR sequences of three Korean PSV strains were in 
the same cluster and were more closely related to the English 
strain V13 than to the Chinese strains csh and YC2011 
(Figure 3A). The 3' UTR sequences of five strains (KS04105, 
KS055217, KS05151, V13, and YC2011) were 82 nucleotides 
long while that of Chinese strain csh was 68 nucleotides in length, 
possibly due to incomplete sequencing (Table S3 in File SI). 
However, they were phylogenetically very close (Figure 3B). 

The leader protein sequences of all PSV strains were 252 
nucleotides (84 amino acids) in length (Table S3 in File SI) and 
show high deduced amino acid identities (95.2-100%) (Table S4 in 
File SI). The PSV leader polypeptide lacked the catalytic residue 
motifs necessary for proteolytic activity and did not contain either 
a zinc-finger motif [Cys2 His2-like fold group] in the leader amino 
terminal region or a tyrosine-phosphorylation motif [KR]-x(2,3)- 
[ED]-x(2,3)-Y]. 

The nucleotide and deduced amino acid sequences of the PSV 
capsid region varied in length from 2430 to 2454 nucleotides 
(encoding 810 to 818 amino acids, Table S3 in File SI). To 
investigate the genetic relationships between the PSV strains, 
pairwise sequence identities were calculated for the deduced 
complete capsid protein sequences of all 6 PSV strains and for the 
sequences of each of the mature capsid proteins, VP 1 - VP4 (Table 
S4 in File SI). All PSV strains showed 88.4% to 97.7% nucleotide 
identities in the complete capsid coding sequences (Table S4 in 
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'(B) 





Figure 1 . Identification and morphology of the porcine sapelovirus (PSV). (A) Immunofluorescence analysis of the PSV infected LLC-PK cells 
by laser confocal microscopy. The LLC-PK1 cells were incubated with mouse anti-153/5B5 monoclonal antibody, followed by staining with 
fluorescein-conjugated goat anti-mouse IgG antibody (green fluorescence). The nuclei were visualized by staining with DAPI (blue fluorescence). (B) 
Electron micrograph (EM) of cultured PSV strain KS05151. Virus pelleted by ultracentrifugation was stained with 1% phosphotungstic acid and 
sprayed onto coated EM grids. 
doi:10.1371/journal.pone.0107860.g001 



File SI). The phylogenetic relationships are shown in Figure 3C. 
The PSV strains had 85.6% to 98.2% nucleotide identities within 
the VP1 coding sequences (Table S4 in File SI). The phylogenetic 
analysis for the VP1 proteins among PSV strains is shown in 
Figure 3D. The VP1 sequences of the KS05151, csh, and YC201 1 
strains encoded an additional 8 amino acids in comparison with 
other strains (Table S3 in File SI). 

The P2 region of all PSV strains was 1989 nucleotides (663 
amino acids) in length and P3 was 2235-2238 nucleotides (745- 
746 amino acids) long (Table S3 in File SI). Pairwise deduced 
amino acid sequence identities of P2 and P3 regions were shown to 
be very high, ranging from 94. 1 % to 99.5% and 95.9% to 99.0%, 
respectively (Table S4 in File SI). 



Analysis of the RNA structures within the 5'UTR 

As described above, the use of 5'RACE allowed determination 
of the complete 5'UTR sequence of the Korean PSVs (some 
491 nt in length); the 5' UTR of these strains was 25 nucleotides 
longer than that of Chinese strain of PSV (YC201 1). This allowed 
prediction of the secondary structure elements within the PSV 
5'UTR (Figure 4). At the extreme 5' terminus were two stem-loop 
structures, labelled domains la and lb. The latter included two 
smaller stem loop structures, labelled Ic and Id. The 5' UTR also 
contained two other domains, labelled domain II and domain III. 
These represent essential components of the IRES and are labelled 
in the same manner as other type IV IRES elements (related to 
that found in hepatitis C virus (HCV), the pestiviruses and certain 
picornaviruses, e.g. porcine teschovirus) [9,30,31]. Domain III 



Table 1. Location of putative cleavage sites in the porcine sapelovirus polyprotein. 





Cleavage site 


Amino acid sequence 


Position of an 


nino acid 








1 a 


2 b 


3 C 


L/VP4 


GNKPQ/GAYNH 


84/85 


84/85 


84/85 


VP4/VP2 


GPSLK/APDKE 


137/138 


137/138 


137/138 


VP2/VP3 


RQ/GFPVR 


375/376 


375/376 


375/376 


VP3/VP1 


YQ/GD 


609/610 


609/610 


609/610 


VP 1/2 A 


AEQL a ' b (ATQT C )/GPYE 


902/903 


894/895 


894/895 


2A/2B 


HDWVQ/GLGQV 


1128/1129 


1120/1121 


1120/1121 


2B/2C 


EPHKQ/GPSDW 


1233/1234 


1225/1226 


1225/1226 


2C/3A 


DAIFQ/GPVQ 


1565/1566 


1557/1558 


1557/1558 


3A/3B 


KQ/GAY 


1665/1666 


1657/1658 


1657/1658 


3B/3C 


KAVVQ/GPDME 


1687/1688 


1679/1680 


1679/1680 


3C/3D 


FVNKQ/GLITE 


1869/1870 


1861/1862 


1861/1862 


3D/ 


F/ 


2331/ 


2323/ 


2322/ 



Letters in bold represent conserved amino acid residues. 
a YC2011, KS05151 and csh strains. 
b KS055217 and KS04105 strains. 
C V13 strain. 

doi:1 0.1 371 /journal.pone.01 07860.t001 
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PI 



P2 



P3 



S'UTR BSMV L VP4 VP2 VP3 
252 159 714 702 



VP1 

855-S79 



2A 2B 
678 315 



2C 3A 3B 3C 

996 300 66 546 



VPgh 



84 53 238 



226 105 



100 22 182 



3D 3"UTR 
1383-1386 



-Poly(A) 



Figure 2. Genome organization of the porcine sapelovirus (PSV). The open reading frames are flanked on either side by UTRs. The numbers 
above or under each rectangle are the length of nucleotides or deduced amino acids. The length of VP1 and 3D regions are different among PSVs. 
doi:1 0.1 371 /journal.pone.01 07860.g002 



contained multiple sub-domains including a pseudoknot (Illf) and 
highly conserved stem-loops Hid and Hie which, by analogy to the 
HCV IRES, interact directly with the 40S ribosomal subunit. 

Analysis of 3'UTR sequence 

The PSV 3'UTR was highly conserved. Based on the size of the 
3'UTR, the PSVs can be divided into a V13-like group (KS04105, 
KS05151, KS055217, YC2011 and V13 strains) which was 82 
nucleotides long and the csh strain that was 62 nucleotides long 
(Figure 5 A). Both of these sapelovirus subgroups possess two 
common domains that can form a stem-loop structure (Domains X 
and Y [Figure 5B-5D]). The V13-like group contained a third 
stem loop region (termed domain Z [Figure 5B, 5C]), located 
upstream of domain Y (Figure 5B, 5C). In the 3 Korean strains 
and YC2011 strain, an interaction of 8 nucleotides between the 
same sequence region within the loop of domain Z (nt 7485 to 
7492 in the KS05151 strain) and the loop of domain X (nt 7549 to 
7556 in the KS05151 strain), forming a loop-loop intramolecular 
"kissing" RNA interaction, appeared to be possible (Figure 5B). 
The VI 3 strain had a potential interaction of 7 nucleotides 
between the loops of domain Z (nt 7411 to 7416) and domain X 
(nt 7475 to 7480), forming a similar intramolecular "kissing" RNA 
interaction (Figure 5C). The csh-like group sequence was not 
predicted to form a loop-loop intramolecular "kissing" RNA 
interaction due to the lack of the domain Z (Figure 5A, 5D). 

Identification of a c/'s-acting RNA element (CRE) 

The CRE is an essential element in picornavirus RNA 
replication [7,32-35]. These relatively short elements can be 
located in different places within the genome; they act as the 
template for the uridylylation of VPg to form VPg-pU-pU and 
contain an essential motif AAAYA [36]. Analyses of the three 
Korean sequences (KS04105, KS05151, KS055217) and other 
known PSV strains including the Chinese (csh and YC2011) and 
English (VI 3) strains revealed the presence of several AAACA 
motifs in the genomes of these strains. As shown in Figure 6, 
however, a stem-loop structure with an exposed loop containing 
the AAACA motif was only found in the 2C coding sequence (i.e., 
nucleotide 4491 to 4550 in the genome of KS05151). This hairpin 
loop included 16 nucleotides. Interestingly, the KS04015, 
KS05105, KS055217 and YC2011 strains had two AAACA 
motifs within this loop structure (Figure 6). 

Discussion 

We report here the isolation of three Korean PSV strains, 
KS05151, KS04105 and KS055217 from porcine diarrhea 
specimens. The Korean PSV strains were identified as PSV by 
RT-PCR, IFA and TEM assays. The genome sequences of the 
Korean PSV strains were determined and proved to be the first 
complete genome sequences for PSVs. They have a genome 
organization typical for members of the genera Cardiovirus , 



Aphthovirus, Erbovirus, Kobuvirus, Teschovirus and Senecavirus 
[1 1]. Moreover, these strains had distinct 2A proteins from those 
of the Enterovirus genus and a 5' UTR with a type IV IRES 
[12,14,15,37]. 

The 5' UTRs of picornaviruses are highly structured and 
contain an IRES that directs RNA translation by internal 
ribosome binding [8,9]. Picornavirus IRES are currently divided 
into five distinct types by the primary sequence, secondary 
structure, location of the initiation codon and activity in different 
cell types [9,10]. In a previous study [16], the IRES elements of 
PSV VI 3 strain and simian sapelovirus SV2 strain were found to 
be related functionally and structurally to the type IV IRES 
element from porcine teschovirus 1 and hepatitis C virus. 
Comparative sequence analysis of the Korean PSV strains with 
PSV VI 3 strain showed that the structural features of the IRES 
elements were well conserved in all PSV species including the 
Korean PSV strains but they lacked a domain IIIc [16]. 

At the 5 ' UTR terminus, enteroviruses and rhinoviruses contain 
a cloverleaf structure which is involved in RNA replication [38] . In 
order to identify whether PSV species have a cloverleaf RNA 
structure at the 5 'UTR, the complete 5 'UTR needs to be known. 
However, the sequences of the known PSV strains including the 
English VI 3, and the Chinese csh and YC201 1 strains lacked the 
5' terminal UU residues, which were necessary for picornavirus 
RNA replication [39]. Using 5' RACE, an additional 25 
nucleotides including 5' terminal UU residues were identified 
compared to the recently sequenced Chinese YC201 1 strain [29]. 
Unlike enteroviruses and rhinoviruses [38], the Korean PSV 
strains had no cloverleaf RNA structure at the 5 'UTR. However, 
two conserved stem-loop motifs were present within the 5'- 
terminal 80 nucleotides (Figure 4). The role of these structures is 
not known but they may be expected to play a role in RNA 
replication analogous to the cloverleaf structure of the enterovi- 
ruses [38,39]. Overall the 5'UTRs of PSVs were quite short, for 
example the poliovirus 5'UTR is about 740 nucleotides in length 
while the foot-and-mouth disease virus (FMDV) 5'UTR is over 
1300 nucleotides [40]. 

The picornaviruses that have a L protein preceding the capsid 
region are members of the genera Cardiovirus, Aphthovirus, 
Erbovirus, Kobuvirus, Teschovirus and Sapelovirus [41]. In 
aphthoviruses and erboviruses, the L proteins are papain-like 
cysteine proteinases that are able to cleave at their own carboxy- 
terminus and also to induce the cleavage of the eukaryotic 
initiation factor (elF) 4G, leading to the shut-off of host-cell protein 
synthesis [42,43]. The L protein of encephalomyocarditis virus (a 
cardiovirus) binds zinc, is phosphorylated during viral infection, 
and has been reported to affect the efficiency of genome 
translation [41]. The properties of the sapelovirus L protein are 
not known; it has neither the catalytic dyad (Cys and His), 
conserved in a papain-like thiol protease found in FMDV L 
protein [44] , nor a putative zinc-binding motif, Cys-His-Cys-Cys, 
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Figure 3. Sequence comparisons and phylogenetic analysis of Korea porcine sapelovirus strains. The phylogenetic tree of 5' 
untranslated region (UTR) sequence (A), 3'UTR sequence (B), P1 nucleotide sequence (C), VP1 nucleotide sequence (D), P2 nucleotide sequence (E), 
and P3 nucleotide sequence (F) were constructed using the neighbor-joining method with 1,000 bootstrap replicates, and the branch length is 
indicated at each branch node. 
doi:1 0.1 371 /journal.pone.01 07860.g003 



found in encephalomyocarditis virus [45]. Further studies are 
required to address the function of PSV L protein. 

In general, picornavirus 3 ' UTRs vary in length between 40 and 
165 nucleotides. The length of the 3' UTR in PSV is 82 
nucleotides (strain VI 3, YC2011, KS04105, KS05151 and 
KS055217) although a shorter sequence (68 nucleotides) was 
described for the csh strain [20]. In a previous report [12], PSV 
strain VI 3 was predicted to include three stem-loop structures in 
the region of the 3'UTR using nine nucleotides of the terminal 
part of the 3D coding region. In the present study, the Korean 
viruses KS04105, KS05151 and KS055217 plus the Chinese 
YC2011 strain each appear to have these 3 stem-loop structures 
(X, Y and Z), but the csh strain showed only two stem-loop 
structures (X and Y) since the sequence was shorter. The domain 
Z, in which the stop codon (UGA) is located, is the most conserved 
region of the 3'UTR within PSVs, whereas domains Y and X were 
considered more variable regions (Figure 5B, 5C), as they show 
heterogeneity in both length and nucleotide sequence. The 
differences between the 5 different PSV strains and the csh strain 
appear attributable to an incomplete 3' terminal sequence for the 
csh 3' UTR region. Moreover, five strains, except for the strain 
csh, were predicted to have an intramolecular kissing RNA 
interaction between the X and Z domains. Due to the lack of 3' 
terminal sequence of csh 3 ' UTR region, no intramolecular kissing 



RNA interaction could be predicted (Fig. 5D). Interestingly, the 
VI 3 strain showed intramolecular interaction of 7 nucleotides, but 
the three Korean and one Chinese strains had intramolecular 
interaction of 8 nucleotides. The 3' UTR plays an important role 
in picornavirus replication. Serial passage of mutant viruses in 
which such interactions were disrupted resulted in production of 
revertants in which the tertiary kissing interaction was restored, 
indicating the functional importance of the interaction in the 
enterovirus 3'UTR [46-50]. Further study is required to 
determine whether these regions are important for PSV replication 
using modifications of the relevant PSV nucleotide sequences. 

A CRE of picornaviruses has been identified in six genera of 
Picornaviridae family; Enterovirus [6,51], Rhinovirus [5], Cardi- 
ovirus [52], Aphthovirus [7], Parechovirus [4], and Hepatovirus 
[53]. However, no putative CRE has yet been reported for viruses 
in the Sapelovirus, Kobuvirus, Erbovirus, Teschovirus, Avihepa- 
tovirus, Senecavirus and Tremovirus genera. The location of the 
CRE in the genomic RNA varies between the picornavirus genera. 
It is located in the coding region for 2C in enteroviruses [6,36], for 
VP1 in species B rhinoviruses [32], for VP2 in species C 
rhinoviruses [38] and cardioviruses [52], for VP0 in parecho- 
viruses [4] and for 3D of hepatoviruses [54]. In FMDV, the CRE 
is located just upstream of the IRES [7,35]. In the present study, a 
putative CRE was located in the 2C coding region of all PSV 
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Figure 4. Sequences and structural features of 5' untranslated region of the porcine sapelovirus KS05151 strain. At the extreme 5' 
terminus are two stem-loop structures, labelled domains la and lb. Secondary stem-loops include two smaller stem loop structures labelled Ic and Id. 
A secondary structure model for the domains II and III of the type IV internal ribosome entry site element is shown (this model is based on previously 
published studies [9,29,30]). 
doi:1 0.1 371 /journal.pone.01 07860.g004 



strains. Generally, the AAACA motif is located in the loop of 
picornavirus CREs and the first and second A residues are 
involved in providing a template for the addition of uridine onto 
VPg [38,39,53,55,56]. In contrast to other known CREs, the PSV 
strains had two copies of the AAACA motif within the sequence 
CAAACATAATAAACAA. This indicated that one or both of 
these AAACA motifs may be involved in being a template for the 



addition of uridine onto VPg. Future functional analyses are 
needed to identify whether one or both motifs are templates for the 
uridylylation of VPg to make VPg-pU-pU 0 H- 

In this study, we characterized the structural features of three 
Korean PSV strains in comparison with the other known PSV 
strains. All PSV strains showed the typical picornavirus genome 
organization. We have identified putative RNA structures in the 



PLOS ONE | www.plosone.org 



8 



September 2014 | Volume 9 | Issue 9 | e1 07860 



Full-Length Genomic Analysis of Korean PSV 



(A) 

csh GAUAUGUUUUAGGCGGCGUGAACAUAUCAGUGAUACAGGAUUAACAAUUAGGCUAAUUGCCAAUAGACCCUAAUCCG 

KS04105 GAUAUGUUUUAGGCGGCGUGAACAUAUCAGUGAUACAGGAUUAAUAAUUAGGCUAAUUGGCAAUAGACCCUAAGCCGCCUAUAGGGUCUAC 

KS05151 GAUAUGUUUUAGGCGGCGUGAACAUAUCAGUGAUACAGGAUUAAUAAUUAGGCUAAUUGGCAAUAGACCCUAAGCCGCCUAUAGGGUCUAC 

KS055217 GAUAUGUUUUAGGCGGCGUGAACAUAUCAGUGAUACAGGAUUAACAAUUAGGCUAAUUGGCAAUAGACCCUAAGCCGCCUAUAGGGUCUAC 

V13 GAUAUGUUUUAGGCGGCGUGAACAUAUCAUUGGCACAGGAUUAAUAAUUAGGCUAAUUGGCAAUAGACCCUAAUCCGCCUAUAGGGUCUAC 

YC2011 GAUAUGUUUUAGGCGGCGUGAACAUAUCAGUGAUACAGGAUUAACAAUUAGGCUAAUUGGCAAUAGACCCUAAGCCGCCUAUAGGGUCUAC 

KSO41O5.KS05151.KSO55217 ^^^^^H I^HHi 



(B) 



7477 
7568 



C \ 

<\\ C 
VA 



• C — - G G ' 
G — c 



(C) 



"g G c 



u *~ 

u / c 



u 



K 



7400' 
7491 



C 
U. 
G^ 



G 



A U 



C 
G. 



U " 

\, G 



>0" 
A 

"C 



Ig c ' 
I 



(D) 



-G G c 



U U 



G 

A Y 

, G G 

7426 A 



7502 



u 



c 

A 



C 
A 
A 



*C-G A 
C-G X 
G — C 



Figure 5. Sequences and structural features of 3' untranslated region (UTR) of the porcine sapeloviruses. (A) The nucleotide sequences 
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5'UTR and 3' UTR plus a CRE in the 2C coding sequence. 
Interestingly, the structural features of the CRE in the 2C coding 
sequence and of the 3 'UTR were different between the strains 
circulating in the recent and past decades. These first complete 
genome data for PSV (for the Korean PSV strains KS05151, 
KS04105 and KS055217) will facilitate future investigations 
concerning the molecular pathogenesis and evolutionary charac- 
teristics of this virus. 

Supporting Information 

Figure SI Phase contrast photomicrographs of control 
and infected LLC-PK1 cells, and RT-PCR assay for 
detecting porcine sapelovirus (PSV) VP1 coding region. 

(A) Mock-inoculated control cells. (B) Cells at 1 day after infection 
with Korean PSV strain KS05 151. Note the shrinking and 
rounding up of the infected cells. Microscope settings Ocular: 10; 
Lens: 10X. Scale bar, 200 |im. (C) RT-PCR with primers specific 
for part of the PSV VP1 coding region generated the expected 
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